Asymptotic convergence of the partial averaging technique 



Cristian Predescu and J. D. Doll 

Department of Chemistry, Brown University, Providence, Rhode Island 02912 

David L. Freeman 

Department of Chemistry, University of Rhode Island, Kingston, Rhode Island 02881 

(Dated: February 2, 2008) 

We study the asymptotic convergence of the partial averaging method, a technique used in con- 
junction with the random series implementation of the Feynman-Kag formula. We prove asymptotic 
bounds valid for most series representations in the case when the potential has first order Sobolev 
derivatives. If the potential has also second order Sobolev derivatives, we prove a sharper theorem 
which gives the exact asymptotic behavior of the density matrices. The results are then specialized 
for the Wiener- Fourier series representation. It is found that the asymptotic behavior is o(l/n'^) if 
the potential has first order Sobolev derivatives. If the potential has second order Sobolev deriva- 
tives, the convergence is shown to be 0(l/n^) and we give the exact expressions for the convergence 
constants of the density matrices. 

PACS numbers: 05.30.-d, 02.70.Ss 

Keywords: random series, Feynman-Ka^ formula, partial averaging, convergence rates 



I. INTRODUCTION 

Path integral methods-'^ are perhaps the most impor- 
tant techniques utilized to account for quantum contri- 
butions in equilibrium statistical simulations for many- 
body systems. They are based on the Feynman-Kag 
formula,'^'- which represents the density matrix and 
hence the partition function of a canonical quantum sys- 
tem as an infinite dimensional stochastic integral of a 
Brownian bridge functional.* This stochastic integral is 
then approximated by a sequence of finite dimensional in- 
tegrals, which in turn are evaluated by Monte Carlo tech- 
niques. The two most important families of discretization 
techniques, while not unrelated, can be distinguished by 
the way the paths are parameterized. The discrete path 
integral (DPI) methods^^Sii utilize the Trotter product 
ruleS and appropriate short-time approximations to con- 
struct a more accurate finite-dimensional approximation 
of the quantum-mechanical density matrix. The random 
series implementations-'^ utilize the Ito-Nisio theorem.^° 
to construct the Brownian bridge in the Feynman-Kag 
formula. Independent of which method is used, a major 
concern in the practical applications is the asymptotic 
rate of convergence as measured against the number of 
variables used to parameterize the paths n as well as 
the computational time necessary to evaluate the vari- 
ous quantities involved. 

The asymptotic behavior of the discrete methods is 
well studiediiiiiSii^ and is known to be 0{l/n^) if the 
expressions of the short-time approximations are func- 
tionals of the potential only. In contrast, the relative 
merits of the random series techniques are poorly un- 
derstood. There are two basic implementations of the 
random series techniques that arc important for practical 
applications: the partial averaging*'^'* and the reweighted 
methodsi^ In this paper, we shall study the asymptotic 
rates of convergence of the partial averaging method for 



different classes of basis entering the Ito-Nisio theorem 
as well as for different classes of potentials. We obtain 
surprisingly powerful convergence theorems and, most of 
the time, we provide not only the convergence order of 
different methods, but also the convergence constants. In 
the remainder of the introduction, we try to convince the 
reader that such an analytical study is important both 
theoretically and practically. 

The partial averaging method has been rarely utilized 
in practical applications because it requires the Gaussian 
transform of the potential for its implementation. For 
real life potentials, this is a difficult but not impossible 
task. However, it was generally considered that the im- 
provement the technique brings in does not warrant the 
effort of computing the Gaussian transform of the poten- 
tial and the so-called gradient corrected partial averaging 
method was used instead. It has been shown that this lat- 
ter method has general 0{l/n^) asymptotic behavior for 
sufficiently smooth potentials and it has been argued that 
there is not much reason to suspect a better convergence 
rate for the full partial averaging methodii^ However, 
more accurate numerical evidence recently presented in 
Ref. iJD suggests that the full partial averaging method 
does have in fact a better behavior: if the technique is 
used in conjunction with the Fourier path integral ap- 
proach (FPI)^*-^^ and if the potential is smooth enough, 
the asymptotic convergence is 0{\/ti?). 

Consistent with these observations, in this paper we 
show that the asymptotic rate of convergence for the 
FPI partial averaging approximation (PA-FPI) is indeed 
0{l/rr') whenever the potential has second order deriva- 
tives in a sense that will be made clear in the text. While 
we shall prove quite general convergence results valid for 
almost all series representations, the PA-FPI approach, 
which is based on the Wiener-Fourier basis, will receive 
a special treatment. This is motivated by the fact that 
the Wiener-Fourier series representation is the optimal 
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representation as to the minimization of the number of 
variables used to parameterize the paths^ 



The importance of the partial averaging method re- 
sides also in the fact that it acts as a prototypical strat- 
egy for improving the asymptotic rate of convergence of 
the random series path integral methods. As such, the 
reweighted random series technique^ achieves superior 
asymptotic convergence by simulating the partial aver- 
aging approach. For this reason, there is a close re- 
lation between the asymptotic rates of convergence of 
the two methods. While the asymptotic behavior of the 
reweighted method will be the object of a separate paper, 
we mention that such a study is greatly simplified once 
the convergence rates of the partial averaging method are 
known. 



The practical importance of the convergence theorems 
that will be established in the present paper consists in 
the fact that they give a priori knowledge about the 
asymptotic order of convergence based only on readily 
verifiable properties of the potential. These theorems 
also provide useful estimates of the absolute error in the 
evaluation of the density matrices and, by integration, of 
the partition functions. Moreover, if the potential has 
second order derivatives, the estimates for the density 
matrix are asymptotically exact and they can be used to 
derive the convergence constants for many other prop- 
erties, as for instance the convergence constants for the 
H-method and T-method energy estimators However, 
such applications require careful analytical work that is 
beyond the scope of the present paper. 



The content of the remainder of this paper is organized 
as follows. In the next section, we give a short review of 
the partial averaging method with the purpose of estab- 
lishing the notation. We also provide a description of 
the main classes of potentials for which the asymptotic 
convergence will be studied. In Section III, we prove the 
main convergence theorems, which are valid for most se- 
ries representations. We prove some asymptotic bounds 
for the density matrix in the case when the potential has 
first order Sobolev derivatives and then we give the exact 
asymptotic behavior for the potentials having second or- 
der Sobolev derivatives. In Section IV we specialize the 
results of Section III for the case of the Wiener-Fourier 
series representation. We prove that if the potential has 
first order Sobolev derivatives, then the convergence is 
o(l/n^). Then, we show that if the potential has sec- 
ond order Sobolev derivatives, then the convergence is 
0{l/n^) and we also provide the exact expression for the 
corresponding convergence constant. We illustrate the 
derived asymptotic convergence constants for the sim- 
ple case of the harmonic oscillator and, in Section V, we 
summarize and discuss our results. 



II. THE PARTIAL AVERAGING STRATEGY 

A. Description of the method 

In this section, we shall give a short review of the par- 
tial averaging method with the sole purpose of estab- 
lishing the notation. For a more complete discussion, the 
reader should consult Ref. (0) and the cited bibliography. 
The starting point is the Feynman-Kag formula 

p{x x'-(3) ^E^^p f'v\xriu)+aB^]du\, 

(1) 

where p(a;, a;'; j3) is the density matrix for a monodimen- 
sional canonical system characterized by the inverse tem- 
perature (3 = l/{kBT) and made up of identical particles 
of mass mo moving in the potential V{x). The stochas- 
tic element that appears in Eq. {B^, u > 0}, is a 
so-called standard Brownian bridge defined as follows: 
if {Bu, u > 0} is a standard Brownian motion starting 
at zero, then the Brownian bridge is the stochastic pro- 
cess {Bu \ i?i = 0, < u < 1} i.e., a Brownian motion 
conditioned on i?i = 0. In this paper, we shall reserve the 
symbol E to denote the expected value (average value) of 
a certain random variable against the underlying proba- 
bility measure of the Brownian bridge B^. To complete 
the description of Eq J^l, we set Xr{u) = x + {x' — x)u 
(called the reference path), a = (/i^/3/mo)^''^, and let 
Pfp{x,x'; f3) denote the density matrix for a similar free 
particle. 

The generalization of the Eq. to a d-dimensional 
system is straightforward. The symbol B^ now denotes 
a d-dimensional standard Brownian bridge, which is a 
vector {B'^ i,B^ 2j ■ ■ ■ d) ^^^^ components being 
independent standard Brownian bridges. The symbol a 
stands for the vector (ai, a2, ■ ■ ■ , <^d) with components de- 
fined by af = ti^(3/mQ^i, while the product aB^ is inter- 
preted as the d-dimensional vector of components cnB'^ ^. 
Finally, x and x' are points in the configuration space 
while Xr{u) ^ X -\- {x' — x)u is a line in connecting 
the points x and x' . In this paper, we shall conduct the 
proofs for monodimensional systems and only state the 
easily obtainable d-dimensional results. 

The most general series representation of the Brownian 
bridge is given by the Ito-Nisio theorem^iS the statement 
of which is reproduced below. Assume given {Afe(r)}fc>i 
a system of functions on the interval [0,1], which, to- 
gether with the constant function \q(t) — 1, makes up 
an orthonormal basis in i^[0, 1]. Let 

Afc(i) = / \k{u)du. 
Jo 

If f2 is the space of infinite sequences a = (ai, 02, . . .) and 

00 

P[a] = II p(ak) (2) 
fe=i 
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is the probability measure on Q such that the coordi- 
nate maps a — > are independent identicaUy distributed 
(i.i.d.) variables with distribution probability 



^(flfc e A) 



1 



then 



= X]afeAfeH, 0<u<l 



(4) 



fe=i 



is equal in distribution to a standard Brownian bridge. 
Moreover, the convergence of the above series is almost 
surely uniform on the interval [0, 1]. 

Using the Ito-Nisio representation of the Brownian 
bridge, the Feynman-Kag formula ^ takes the form 

= / dP[a]exp|-/3 /Vk(.) 
Pfp[x,x';(J) Jn L Jo L 



CT^afeAfe(u) 



(5) 



The independence of the coordinates Ofc , which physically 
amounts to choosing those representations in which the 
kinetic energy operator is diagonal, is the key to the use 
of the partial averaging method. Denoting by E„ the av- 
erage over the coefHcients beyond the rank n, the partial 
averaging formula reads 

X exp 



d^(ai) . . . / d^(a„) 



13 E„ / V 





.{u) + a'^akAkiu) du>. (6) 



fc=i 



To make more sense of the above formula, it is conve- 
nient to use the notation introduced in Ref. Q) 

n oo 

'S'"(a) = ^afcAfc(M) and B'^ia) ^ ^ afeAfc(u), 



k=l 



k=n+l 



which denote the nth order partial sum and and nth or- 
der tail series, respectively. By construction, the random 
variables 5" (a) and i?"(a) are independent and 

oo 

fc=i 

Moreover, the variables and have a joint Gaussian 
distribution of covariances 

'DC OO 

A;— n+l k—n-\-l 

oo 

and 7„(u,t) =E„(B:b;?) - ^ Afc(w)Afc(T). (7) 

k=n+l 



Equivalently, by using the fact that S*" (a) and B" (a) are 
independent and their sum is B'^(a), one may evaluate 
the above series to be 



k=l 



(3) En{B:r = E{Bir - ^ A^iur = u{l - ^) - ^ Ak{u)\ 

n n 
k=l 

and 

n 
k=l 

n 

= min(u, r) - ur - E Afe(u)Afc(r). 



fc=i 

The computation of E{B°f, HB°)^, and E{B^B°) is 
straightforward and is performed in Appendix A. 

Going back to Eq. 0, one inverts the order of inte- 
gration in the exponent and computes 

E„ / duV[xr{u) + <7Bl{a)] 
Ja 

= f duEnV[xr{u)+aS'^{a)+aBll{a)] (8) 



where 



exp 



2r^(«) 



V{y + z)dz, 



with r^j(u) defined by 



Tliu)=a^E,,{B:f^a' 



(9) 



fe=i 



(10) 

In deducing Eq. ©, one uses the fact that the variable 
has a Gaussian distribution of variance E„(i3^)^. To 
summarize, we define the nth order partial averaging ap- 
proximation to the diagonal density matrix by the for- 
mula 



Pfp{x,x';(3) 



d^(ai) . . . / d^(a„) 



X exp 



13 / Vu,. 

'0 



^(u) + crE«feAfe(u) dwWll) 



k=l 



It has been shown- that the rate of convergence of the 
partial averaging method is controlled to a first approx- 
imation by the rate of decay of T'fi{u) to zero. More 
precisely, the value T^{u)du may be taken as a mea- 
sure of the efficiency of different basis sets {Afc(T)}fc>i. It 



4 



was found that T'^{u)du attains its unique minimum understand 



on the cosine-Fourier basis {Afc(T) — -y/2 cos(/c7rT); k > 
1}. The corresponding series representation for the 
Brownian bridge is the Wiener- Fourier construction 



2 Z-^ 

k=l 



a-k- 



sin(fc7ru) 
1 



< u < 1 



and was used by DoU and Freeman*^ in their definition 
of the Fourier Path Integral method. Though their work 
was based on arguments other than those presented here, 
we consider that the name of Wiener-Fourier Path Inte- 
gral (WFPI) method is more appropriate because the 
random series 



aou 



I 2 sinffcTTu) 



k=l 



k 



was historically the first explicit construction of a stan- 
dard Brownian motion and was due to Wiener^ 



B. Mathematical considerations 

In this subsection, we discuss the classes of basis 
{Afc(u)}fe>i as well as the classes of smooth enough po- 
tentials for which we study the asymptotic convergence. 
Anticipating later results from this paper, the rate of 
convergence of the partial averaging sequence of approx- 
imations of the density matrix will be shown to depend 
solely on the properties of the function 7„(it, r) defined 
by Eq. Q. Regarded as an integral kernel, the func- 
tion 7„(w, r) defines a positive definite quadratic form 
on L'^([Q, 1]), that is 



7„(w, T)f{u)f{T)AuAT > 0, V/ e L2([0, 1]). 



^0 



It is but a simple exercise to verify that in general 
7„(m, r)'^ defines such a positive definite quadratic form 
for all integers fe > 1. In particular, we are interested in 
the convergence properties of 7„(w, r)^. 

In this paper, we shall restrict our discussion to those 
basis {\k{u)\k>i for which there is a distribution g2{u, r) 
such that (i) 



lim 



/oVoNn(M,T)2dud7 



(12) 



in the sense of distributions and (ii) g2{u,T) defines an 
integral kernel on ^^([0,1]) that is strictly positive def- 
inite. By convergence in the sense of distributions, we 



lim 



92{u, T)h{u, T)dudr 



"'0 



for all continuous functions h(u, r), while the strictly pos- 
itive definiteness of g2{u,T) is the statement 



92iu,T)f{u)f{T)dudT = 



/ = a.s. 



We mention that it is not clear if in fact these two 
conditions are satisfied for all basis {Xk{u)}k>i- How- 
ever, the reader may always verify them for the problem 
at hand. We prove in Appendix B that for the Wiener- 
Fourier basis 



lim 



/oVo^7n(u,T)2dudr 



S{u — t) 



(13) 



in the sense of distributions. Moreover, 6{u — t) is indeed 
strictly positive definite because 



S{u ~ T)f{u)f{T)dudT 



fiu)^du 



is zero if and only if / = a.s. 

A special category of basis {Xkiu)}k>i are those for 
which the primitives Afc(u) do not change sign. An 
example is furnished by the Haar wavelet basis, which 
generates the so called Levy-Ciesielski representation of 
the Brownian bridge>iii^ For such basis, the functions 
7„(m,t) are positive and we shall further assume that 
there is a distribution gi{u, r) such that (i) 



lim 



(14) 



in the sense of distributions and (ii) gi (u, r) defines an 
integral kernel on ^^([0, 1]) that is strictly positive defi- 
nite. We analyze this category of basis separately because 
there are more general convergence results available for 
them. 

In the second part of this subsection, we shall dis- 
cuss the class of potentials for which the convergence 
is analyzed. For sure, the analysis must be restricted 
to those potentials for which the Feynman-Kag formula 
holds and for which the method does in fact converge. 
The Feynman-Kac formula Q is known to hold for a 
quite large class of potentials which is called the Kato 
class. On the other hand, it has been proven^^ that the 
partial averaging method is convergent for all series rep- 
resentations and for all Kato-class potentials that have 
finite Gaussian transform. We provide below the mathe- 
matical definition of this class. The readers who find this 
definition meaningless might get a better insight about 
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the nature of the potentials in this class by studying the 
examples we later provide. To arrive at the definition of 
the Kato class^ii we let 

r \y\ d = i, 
9{y)^\ H\\y\\-') d^2, 

and define the Kato class Kd as the set of all measurable 
functions / : R'* ^ M such that 

lim sup / \f{y)g{x - y)\dy = 0. (15) 

"iO xeS.'' J\\x-y\\<a 

We also say that / is in K^^"" if Id/ e Kd for all bounded 
domains D C M''. We say that V{x) is a Kato-class 
potential if its negative part = max{0, ~V} is in Kd 
while its positive part V+ — max{0, V} is in K]^^. 

Roughly speaking, Eq. H15|l imposes a restriction on 
the singularities of the potentials. As far as the chemical 
physicist is concerned, the standard prototype as to the 
nature of the singularities is the general n-body potential 

l/(r-i,...,r-„)=.f]— (16) 

Z— 1 " " J 

where the individual particles are assumed to move in 
the usual tridimensional space. Using the observation 
that both Kd and K^^'^ classes are linear spaces and also 
using Proposition 4.3 of Ref. the reader may prove 
that the potentials of the form given by Eq. H16() are 
in the class Kd for all g < 2. Such potentials include 
for instance the coulombic potentials as they appear in 
electronic structure calculations and for which q = 1. 
Remark that the singularities in Eq. H16() can be oriented 
upward or downward. It does not make any difference 
as to the validity of the Feynman-Kag formula as well as 
the convergence of the partial averaging method. 

Another important example of Kato-class potentials 
are the functions V{x) that are continuous and bounded 
from below. It is trivial to verify that if ci > and 
C2 > are some nonnegative numbers and Vi{x) and 
V2{x) are Kato-class potentials, then ciVi{x) + C2V2{x) 
is also a Kato-class potential. Using this observation, one 
may argue that any ab initio potential obtained at the 
level of the Born-Oppenheimer approximation must be 
a Kato-class potential. This is so because it can usually 
be written as a sum between a Kd potential of the type 
given by Eq. ()16fl (the internuclear repulsion term) and a 
continuous and bounded from below potential (this term 
arises from the electron-nuclear and electron-electron in- 
teractions). More generally, even for strongly ionic sys- 
tems, one notices that the Born-Oppenheimer potential 
can have at most coulombic negative singularities and 
therefore it is of Kato-class. 

Many of the properties of the Kato-class potentials 
were enumerated by Aizenman and Simoni^ As argued 



by them, the Kato-class is the natural and in many as- 
pects maximal class of potentials for which the Feynman- 
Kag formula holds. It is known for instance that any 
Kato-class potential is locally integrable. That is, 

/ \V{x)\dx^ I lr,|T^(a;)|da; < cx) 

for all bounded domains D G R'^. Because of this re- 
striction, the Leonard- Jones potential as well as some 
other empirical potentials are not included in the Kato- 
class due to the r~^^ singularity. Even if the Feynman- 
Kag formula may be well defined for such potentials with 
positive non-integrable singularities, it is not clear if the 
computed answer is the Green's function of the corre- 
sponding Bloch equation. However, as far as the chem- 
ical physicist is concerned such issues are hardly of any 
relevance because the Kato class comprises all potentials 
of physical interest. The Leonard- Jones potential can be 
brought into the Kato-class by truncation or other ap- 
proximations. 

The condition that the potentials have finite Gaussian 
transform can be formulated as follows. For all vectors 
a — («!,..., Qfd) £ M.'^ of strictly positive components, 
we consider the Gaussian measure 

exp(^-^^jdzi...dz.. 

Then, a potential is said to have finite Gaussian trans- 
form if 

Ma(a;)=/ \y{x + z)\Aii^{z)<^, (18) 

for all a; e M'' and a G M'j^. Let us notice that from the 
practical point of view, this additional condition is hardly 
a restriction. The potentials used in actual simulations 
are usually a sum between the Born-Oppenheimer inter- 
atomic potential and a constraining potential. The first 
one always has a finite Gaussian transform. The second 
one is added in order to ensure that the partition function 
is finite 

Z(/3)<oo, V/3>0 

and is intended to simulate, for example, the effect of 
the container in which a reaction takes place. For this 
purpose one usually uses a continuous and bounded from 
below potential, which is, of course, a Kato-class poten- 
tial. If, for instance, the potential increases to infinity 
slower than any exponential exp(c||a;|p) with c > 0, then 
it has a finite Gaussian transform. But this is hardly a 
restriction since the typical constraining potentials are 
polynomials. 

Under these conditions, it was showniSi that the partial 
averaging sequence of approximations is convergent to 
the correct Feynman-Kag result (which by itself is prop- 
erly defined) at least as far as the convergence of the 
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density matrix and the partition function is concerned. 
However, we expect that the rate of convergence depends 
on the smoothness of the potential. As mentioned in the 
introduction, we are interested in establishing the class of 
potentials for which the fastest convergence is achieved. 
While the maximal class depends upon the specific series 
representation, we shall see that there is a sufficiently 
large class of potentials for which all series representa- 
tions achieve their fastest asymptotic convergence. 

Theorem 3.b) of Ref. (jH)' says that in order to verify 
that a potential has a finite Gaussian transform, it is 
enough to verify that 



I^L(^o) 



\V{xo + z)\diJ.a{z) < oo, 



(19) 



for all a G R^f. and an arbitrarily given xq G M.'^. If we 
choose Xq = 0, then the set of potentials having finite 
Gaussian transform is the set naI/^(M'^), where Ll^{M.'^) 
is the space of functions / for which the weighted norm 



l/ll 



\fiz)\d^,M) 



(20) 



is finite. By Theorem 3. a) of Ref. (HI), ^^(M'*) C 
Lj^^CR''-) and conversely, one may argue that a function 
has finite Gaussian transform if it is locally integrable 
and its modulus increases slower than any Gaussian at 
infinity. 

However, in this paper we shall assume that the po- 
tential V{x) lies in the set n^L^ (R'^), where L^(R'^) is 
the space of functions / for which the weighted norm 



\f\[ 



1/2 



(21) 



is finite. Since the measure dfia{z) is a probability mea- 
sure, one may apply the Cauchy-Schwartz inequality and 
see that < ||/||2.a, so that Ll{R'') c ^^(R'') 

and n^i^^Rd) ^ r\aLl{R'^). Therefore, the class of 
potentials discussed in the present paragraph have fi- 
nite Gaussian transform. We remind the reader that 
Ll(R'^) C Lf^M'^) and that the spaces ^^(R'^) are 
Hilbert spaces. 

Still, the class of potentials introduced in the previous 
paragraph is not smooth enough for the purpose of study- 
ing the asymptotic convergence. Later in the paper, we 
shall see that natural classes of smooth enough potentials 
are intersections of weig hted Sobolev spaces n^W^^^iR'^) 
for m = 1 or TO = 2. Since the functions f{x) G L^(R'') 
are locally integrable, they have partial distributional 
derivatives of any order. The weighted Sobolev space 
W™'2(R'') is the space of all functions / e L^(R'') whose 
partial distributional derivatives up to the order to are 
also L^(R'') functions. In an equivalent definition, one 
may argue that / G r\aW^-^^{R'^) if and only if 




/(x)exp ->^-^ GW^ 



for all a G R^, where W'^''^(R'^) is the usual Sobolev 
space on R"^. 

There are two main categories of Kato-class poten- 
tials which we study in this paper: those that also lie 
in HaW^''^ {R''') and, more restrictively, those that lie 
in nQW^^^(R'^). For the first class of potentials, the 
weighted Sobolev norms 

J^^^vixf +J2[d,vix)r^d^,^ix) 

are finite, while for the second class of potentials, the 
norms 

/ I vixy + J2i9.vix)r + J2id^,,vix)r \ d^^^ix) 

are finite. These requirements imposes further restric- 
tions on the singularities of the Kato-class potentials. We 
leave it for the reader to show by explicit computation 
that the potentials of the form given by Eq. H16|l lie in 
the naWyiR^"-) space if and only if g < 1/2 and in 
the naW^'^(R^") space if and only if q < -1/2. For an- 
other example, a general 3n-body potential of the form 
given by the Eq. Hlt)|l but where the inter-particle poten- 
tial is replaced with the Morse potential lies also in the 
naM^^'^(R^") space. 



III. ASYMPTOTIC CONVERGENCE OF THE 
PARTIAL AVERAGING METHOD 

In this section, we shall establish several general con- 
vergence theorems valid for all series representations and 
for smooth potentials. While for most of the paper we 
adopt a rigorous mathematical style, we shall not pro- 
vide a formal proof of the key conjecture expressed by 
Eqs. (|26|l and H27|) . Instead, we shall try to present argu- 
ments to support this plausible assertion within reason- 
able doubt and then explore its implications. 

We start by defining 



Un{x,x' ,f3;a) ^ j Vu,n[^r{u) -I- Q-y^ QfcAfc(u)]dM 
•^0 k=l 



and 



„1 oo 

Uoo{x,x' , f3;d) = / V[xr{u) + cry2akAk{u)]du. 

The variables x, x' , and (3 are interpreted here as pa- 
rameters and, just as a reminder, we separate them from 
the "true" variable a by a semicolon. We notice that 
Un{x,x',(3;a) as a function of a depends only on the 
first n coefficients Oi, . . . , a„. Moreover, in terms of these 
functions, the relation (jHJ reads 



Un{x,x' , (3;a) = E„ [Uac{x, x' , (3; a)] 



(22) 
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and in fact, 

lim Un{x,x',P;a) ^Uooix,x',P;a), (23) 

n — *oo 

as shown by Theorem 1 of Ref. |T^ . To continue with 
the introduction of the notations, we define 

Xnix,x',P;a) = p/p(a;,x';/3)exp[-/3[/„(a;,x',/3;a)] 

and 

Xoaix,x' , (3;a) = pfp{x,x';P)exp[-pUoo{x,x',(3;a)]. 
Then, we have 

p{x,x';P)^E[X^{x,x',f3;a)] (24) 

and 

pPA(x,x';/3)=E[X„(a:,a:',/3;a)], (25) 

respectively. 

A little algebra shows that 

E„Xoo(a;,x',/3; a) - Xn{x,x' ,(3; a) = Xn{x,x' ,(3; a) 

However, for large n, we can expand the exponential in a 
Taylor series and retain the first non-vanishing positive 
term in the series, which is also the one that controls the 
asymptotic convergence. We have: 

E„Xoo(a;,a;',/3;a) - Xn{x,x' , I3;a) « X„(x, x', /3; a) 

X y E„ [U^{x, x',13- a) - Unix, x' , P; a)]' . (26) 

From now on, if An > and i?„ > 0, we interpret An ~ 
Bn to mean 

lim An/Bn = I, 

n — >oo 

while An ~ Bn is interpreted to mean 
limsupA„/_B„ < 1. 

n — >oo 

The error in the equation ]2<i\i is of the order 



3! 



En \Uooix,x' , P; a) - Un{x,x',P; a)]" 



which decays at a faster rate than the local variance 
E„ [Uoa{x, x', P; d) - Unix, x' ,P; a)]^ . 

Therefore, the use of the symbol ~ is justified any time 
the local variance is not zero. However, this property is 
true for almost all a whenever the potential is not con- 
stant. In Eq. H26|l . the term of order one in the Taylor 
expansion cancels because of the identity H22|l . so the 
asymptotic behavior is dictated by the local variance of 



the function Uoa{x,x' , P;d). The Eq. is expected 
to be true for all potentials V{x) S nc<i^(M) which are 
not constant and, from now on, we shall assume that the 
potential satisfies these criteria. Taking the total expec- 
tation E in Eq. (|26|l and remembering that Xn{x, x' , P; d) 
does not depend upon the coefficients beyond the rank 
n, we obtain 

p{x, x'; P) - p^,^{x, x; P) k. — E {X„(x, x' , P; d) 

X \U^(x,x' ,P;d)-Un{x,x' ,P;d)f). (27) 



A. Potentials having first order Sobolev derivatives 

If V{x) € r\aWy^), Eq. inn and Eq. of Ap- 

pendix C show that 

E„r,'(a;,x',/3; a) < En\U oo{x, x ,P;d) -Un{x,x ,P;d)f 
< EnTnix,x',P;d), (28) 

where T„(x, x' , P; d) and Tn{x, x' , /?; a) are positive func- 
tions defined by the equations 



Tn{x,x',P;a) ^ a / du dr jn{u,T)V^ n[^r(u) 
Jo Jo 

+aS:id)]v'^^n[Mr)+aS-{d)] (29) 

and 



Tn{x,x' , P;a) = a I du j dr ^n{u,T)V'[xr{u) 
Jo Jo 

+aBl{d)]V'[xr{T)+<jB'l{d)], (30) 

respectively. Here, V'{x) denotes the first order deriva- 
tive of V{x), while „(x) is the first order derivative of 

Vu,n{x). 

With the help of the functions T,' (x, x', /3; a) and 
Tn{x, x' , P;d) and of the equation (|26|l one may write 



P' 



E„ [Xnix,x',P; a)T^{x,x',P;a)] 



^ E,nX^{x,x',P;a) - Xnix,x' , P;a) (31) 
p2 

^ -^E„ [Xn{x, x' , P; a)Tn{x, x' , P; a)] , 

where we used the fact that the function Xn{x,x' , P;d) 
does not depend upon the coefficients beyond the rank 
n and thus, it can be placed inside the En sign. Now, 
taking the total expectation E in the above equation, we 
obtain 



/32 

[Xnix, x', P; a)Tn{x, x',P; a)] 

<p{x,x';P)-pl^{x,x';P) 
P^ 

< ^E [Xn{x, x',P; d)Tn{x, x', P; a)] 



(32) 
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Moreover, due to the relation (|23|l . we have 

AX„(x,x',/3; a) = X„(x, a;', ^; a) - Xoo(a;,a;', ^; a) 

as n — > cxD and we expect that 

^.^ E [AXn(a;, a:^ /3; d)T„(x, x' , P; a)] _ ^ 
n^oo E[X^{x,x' , P;a)Tn{x,x' , /3;a)] 

because the function Tn{x , x' , (3] a) is positive while 
Xoo{x,x' , (3]a) is strictly positive. Therefore, the decay 
of AX„(x, x' , P; a) to zero makes the numerator converge 
to zero at a faster rate than the denominator. A re- 
lation similar to the one above holds for the function 
T^{x, x' , (3; d) and consequently, we may replace the esti- 
mate (|32ll with 

/3 



where 



[Xoo{x, x', P; a)T^^{x, x' , P; a)] 
<p{x,x';p)-pl^{x,x';p) 
< —E[X^{x,x',P;a)T„{x,x',P;a)] , 



(33) 



With these preparations, we are able to state the fol- 
lowing theorem: 

Theorem 1 Assume V{x) is a Kato-class potential that 
lies in HaW^iR). Then, 

p{x,x';P)- p^^ix,x';p) 

<^ / du/ dr7„(u,r)<,,(«,r), (34) 



K' 



[^,iu, r) = a^E^X^ix, x' , P; d)V'[xr{u) + aBl{d)] 

xy'[a;,.(r) +aB°(a)]}. (35) 



If in addition 7„(u,r) > 0, then the following stronger 
result holds 



p{x,x';P) - p„ {x,x';P) 



P^ \ f f 

— du dT7„(w, r) 

^ L"'o Jo 

1 ^1 

du dTgi{u,T)K^ ^,{u,t) 
Jo ' . 



(36) 



Proof. The first part of the theorem follows directly 
from Eqs. (EH and by mere substitution. In order 
to prove the second part of the theorem, one divides the 
terms in Eq. H33|) by 



1 .1 

du dTJn{u,T) 

Jo 



and then uses Eq. 1)14(1 as well as the expressions 1(29(1 and 
((30(1 to show that 



P^ E[X^ix,x',P;d)T^ix,x',P;d)] 

— lim 

2 ™ /odu/odT7„(«,r) 

P^ E[Xooix,x',P;d)Tr,{x,x',p;d)] 

— lim — 

2 n — >oo 



= lim 

n — ^oo 



/odu/o^dr7„(u,T) 



p{x,x';p) - pI^{x,x';P) 
J^duJ^dTj„{u,T) 



du / dTgi{u,T)Kf^^ 
Jo 



(37) 



In this respect, one should notice that the use of Eq. 1(14(1 
is justified by the fact that x'i'^^^) ^ continu- 
ous function in the arguments (m,t). By construction, 
^,{u,t) is symmetric under the change of variables u 
and T, so it is enough to restrict its analysis to the case 
T < u. However, if r < u, we have 

Jr Jr 

xp[z, x'; (1 - u)P]V'{y)V'{z) dydz. (38) 

The continuity of ^,{u,t) then follows from the con- 
tinuity of the density matrix p{x,x']P) in all arguments 
and from the fact that the integrals appearing in Eq. I(38() 
are finite if ^(a;) e naW^iR). 



Equation 1(37(1 implies the relation 1(36(1 provided that 
the common limit is not zero or infinity. We have 



[ du [ dTgi{u,T)K^^, 
Jo Jo 



(w, r) < supiff ^,(-u, r) < oo, 



where we used the fact that the integral of gi(u, r) on the 
set [0,1] X [0,1] is 1 as well as the fact that K^^,{u,t) 
is continuous on the same (compact) set, thus bounded. 
On the other hand, using Eq. I(35() and the fact that the 
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kernel gi{u,T) is strictly positive definite, the equality 

X f du f dTgi{u,T)V'[xr{u) + aB°{a)] 
Jo Jo 

xV'[xr{T)+aB^{a)] \ = 



implies V'[xr{u) + aB^{a)] = for almost all a. However, 
this is not possible because it contradicts the fact that 
the potential V{x) is not constant. Therefore, the limit 
in Eq. H37|) is finite and non-zero and the proof of the 
theorem is concluded. □ 

Let us notice that the asymptotic rate of convergence 
for the partial averaging method depends upon the na- 
ture of the basis {Xk{u)}k>i solely through the behavior 
of the kernels 7„(u,t). The kernel K^^,{u,t) is inde- 
pendent of the particular basis and has the important 
property 



(39) 



The relation ^ can be proven as follows (without loss 
of generality, one may assume t < u): 

-u,\~t)^cj^ I I p[x,y- (1 - u)P] 

xp[y, z- [u - r)/3]p(z, x' ■ r/3)y (y)y(i) (z) dy dz 

p{x', z; T/3)p[z, y; {u - T)/3]p[y, x; (1 - u)f3] 

xT/«(y)t/(i)(z)d2/dz = K^,Jt,u) = K^,,Ju,t), 

where we used the symmetry of the various density ma- 
trices appearing in the above formula. 

The relation H35|l remains true for d-dimensional sys- 
tems provided that we set af = It? (3/ (2mo,i) and redefine 
the kernel iff ^, (u, r) to be 



p{x,y]T(3)p[y,z; {u - 



xp[z,x';il-u)P] i^Y^^aUd^Viy)d,Viz)]^ dydz.m 



B. Potentials having second order Sobolev 
derivatives 



Theorem 1 provides a useful bound for the asymptotic 
rate of convergence, but gives the exact asymptotic con- 
vergence only if 7„(u, r) > 0. In general, if these last 
functions are not positive, additional cancellations might 
occur if the potential is smooth enough. In this case, 
we need to establish lower and upper bounds that are 
sharper than those offered by Eq. 



It turns out that the natural class of potentials on 
which all series representations attain their fastest con- 
vergence is the class naVFj'^(K). For potentials V{x) in 
this class, Eqs. (|C8ll and IjClOp of Appendix C provide 
the following bounds 

EnY^{x,x' , (3;a) £ En[Uaoix,x' , P;d) - Un{x,x' , P;d)f 



< E„r^'(x,x',/?;a), 



(41) 



where Y^{x, x' , (3] a) and Ynix, x' , a) are positive func- 
tions defined by the equations 



4 /•! A 



Y,^{x,x ,(3;a) ^T„{x,x ,P;a) + — du dT7„(u,T) 



"'0 



xV^^Xriu) + aS:{d)]V,Jxr{T) + aS:^{d)] (42) 



and 



y,:'(x,x',/3; d) = -K(x,a;',/?;a) + r„(x,a:',/?;a)], 

(43) 

respectively. Here, Vi „^{x) denotes the second derivative 
of the averaged potential Vt,n{x). By replacing the func- 
tions Yj^{x , x' , (3; d) and YJ/ {x , x' , (3; d) in Eq. and 
taking the total expectation E, one deduces the following 
sharper version of Eq. (|32|l 



E[Xn{x,x',(3;d)Y;{x,x',f3; d)] 



<p{x,x';(3)-pI^{x,x';P) 
<^E[X„(a:,a:',/3;a)r^'(x,a;',/3;a)] 



(44) 



In deducing Eq. (|44|l . one should remember that the 
functions Xn{x,x' , P;d) do not depend upon the coef- 
ficients ftfe with k > n + 1 and, therefore, E [X„y^] = 
E{X„E„K]} and E [XnY^'] = E{X„E„K']}. 

In the next few paragraphs we prove that the signs ^ 
in Eq. H44|l can be replaced by ~ and then we prove an 
analog of Eq. (jSHl- As shown by Eqs. (|C11|I and l)C12(l 
of Appendix C, 

^4 

EnY^'{x,x' , (3;d) -EnY^{x,x , (3;d) < — du dr 



X 7„ (w, t)2 I E„ { V" [xr (u) + a El {d)]V" [x^ (r) (45) 

Using this relation together with the Eq. (|12|l , one readily 
proves that 



lim 

n — *cxj 



E[X„(x,a;^/?; d)Yll[x,x' ,(3- d)\ 
/odw/pdr7„(u,T)2 



E[X„(x,x',^; a)y^(a:,a:',/3; a)] 
/o^du/ndT7„(w,r)2 
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The last relation implies 

E[X„(a;,x',/3;a)r^'(a;,a;',/3;a)] 

« E [X„(x, x\ (3- a)Y:^{x, x' , /?; a)] , (46) 

provided that 

E[X„(x,x^/j;a)y^(x,x^/3;a)] ^ ^ ^^^^ 
/gdu/Qdr7„(u,T)2 

In general, if a„ > 6„ > are some sequences of real 
numbers, such that 

hm (a„ - 5„) =0 

n — >oo 

and liminf„ 6„ > 0, then a„ sa 6„. Indeed, 

, ^ a„ On - fon ^ , , lhii„(a„ - &„) 

1 < lim — = 1+ hm — < IH — — = 1. 

rwoo 5„ n^oo 6„ hm ml On 

Now, to prove the relation (|47|l and therefore Eq. (|46|l . 
one uses Eq. (|12|l together with Eq. (|42|l to show that the 
expression (I47|l is greater or equal than 

^E|xoo(a;,a;',/3;a) J ^'^J '^t92{u,t) 
xV"[xr{u) + aB°{a)]V"[xr{T) + aB^{d)] 

Because 32(^5 t) was assumed to be strictly positive defi- 
nite, the last quantity is zero only if V"[xr{u)+aB^i^{a)] = 
for almost all a, which only happens when the potential 
V{x) is linear (i.e., its second order derivative is zero). 
Therefore, the relation H46|l holds for all non-linear po- 
tentials. However, the relation H46|l holds for linear but 
non-constant potentials too, because in this case 

EnY^'{x,x' , (3;d) = EnY^{x,x' , f3;d), 

as shown be Eq. H45|) . Therefore, the signs ^ in Eq. H44|) 
can be replaced by «. 

We now use the positivity of the functions 
Y^{x, x' , /3; a) and YJ^'{x, x' , (3; d) and argue as in the proof 
of Eq. (|33|l to conclude that 

E [Xn{x, x', P; d)Yj^{x, x' , P; a)] 
^ E [Xoo{x,x' , (3; d)Y^{x,x' , (3; a)] 

and 

E[Xn{x,x',P;d)Y:^'ix,x',P;d)] 
« E [Xooix, x\ 13; d)Y^'{x, x' , (3; d)] , 

respectively. Corroborating the last two relations with 
Eq. and Eq. (gH), we learn that 



^E[X^{x,x',(3;d)Y:^{x,x',(3;d)] 
^ p{x,x';P)- pI^{x,x';P) 
a — E[Xoo{x, x\ P; a)Y" {x, x' , l3; a)] . 



(48) 



With the help of Eq. H48|l . we are ready to prove the 
following theorem 



Theorem 2 Assume V{x) is a Kato- class potential that 
lies in r\aW^'^{R). Then, 



p{x,x';P) - p^'^{x,x';P) ^ '^'^ J dT7„(M,r) 



x<.-Kr)-- 



1 .1 

du / dT7„(u, r)^ 
"'0 
1 ri 

du / drg2 {u, t)QI ^, {u, r) 
Jo 



(49) 



where 



^,{u, t) = a*E{Xoo{x, x', P; d)V"[xr{u) + aB°(a)] 

xV"[xr{T)+<jB°{d)]y (50) 



Proof Let us notice that if a„ w &„, then a„ w 26„ — a„ 
Indeed, 



hm — 2 iim — 



1 = 1. 



Using this observation together with the equation 148|l . 
we deduce that 

pix,x';P)- pf,^ix,x';P) 

« ^E [Xoo{x, x', p; d)Ynix, x',P; a)] , (51) 

where 

Yn{x,x',P;d) = 2Y^'{x,x',P;d) - Y^{x,x' , P;d). 
From Eq. and Eq. we learn that 

Ynix,x\P;d) = Tn{x,x',P;d) - Z„(a;, a;', /?; a), 
where 

^4 <.! ^1 
Zn{x,x' ,P;d) = — du dr7„(u,T)^ 
^ JO Jo 

xVljxriu) + aS:id)]v'^Jxr{T) + aS^id)] 

If we also introduce the random variable 

^4 



Z'n{x,x',P; a) 



1 

du / dr7„(u,T)^ 
Jo 



du dTg2{u,T)V"[Xr{u) + (jBl{d)] 

Jo 



•KV"[Xr{T)+aB°{d)] 



we notice that the right-hand side of Eq. 149|) is precisely 



P^ 



E{X^{x,x' ,P;a) [Tn{x,x' ,P;a) - Z'^{x,x\P] a)]} . 



Now, if the potential is linear but not constant, the the- 
orem follows trivially because the functions Zn{x,x' , P]d) 
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and Z'^{x, x' , /?; a) are identically zero. Therefore, for the 
remainder of the proof, we assume that the potential is 
not linear. Let us set 

AZn{x,x' , (3;a) = Zn{x, x' , P; d) - Z'^{x,x\ (3]d) 

and notice again that the theorem follows trivially from 
Eq. H51(l provided that 

|E [Xoa{x, x',P; d)AZn(x, x' , (3; d)]\ ^ ^ ^^^^ 
oo E[Xocix,x' , P;d)Yn{x,x\ P;d)] 

In order to prove that the limit in the last equation is 
indeed zero, we observe that since E [Xooi^n] ~ E [Xqo^^] 
and since y,' > Z„, we have 

^.^ E [Xoo{x, x\ /3- d)Yn{x, x\ l3; d)] ^ ^ 
n^oo E[Xoo{x,x' , (3;d)Zn{x,x' , /3;d)] ~ 

Therefore, the limit in Eq. (|52|l is smaller or equal than 

|E [X^{x, x',(3; d)AZn{x, x',(3; d)]\ 
n->oo E[Xoc{x,x',P;d)Zn{x,x',(3;d)] 

However, the limit in Eq. (|53|l is zero because 
E [X^Zn] ~ E [X^Z'„]. This last statement follows from 
the equality 

E [Xooix, x',(3; d)Zn{x, x' , /3; a)] 



lim 

n — ^oo 



lim 

n — *oo 



/„Mu/gdr7„(u,T)2 
E [Xocix, x', 13] d)Z'^{x, x' , (3; d)] 
/gdu/gdr7„(u,T)2 

= du^ drg2(w,r)gf_^,(M,r), 

because the above common limit is finite and non-zero 
by an argument similar to the one employed in the proof 
of Theorem 1. The proof of Theorem 2 is therefore con- 
cluded. □ 

The kernel ^,(u,t) shares most of the properties 

proved for ^,{u,t). In particular, Q^j.xi{utt) is sym- 
metric at the permutation of the variables u and r and, 
if T < w, 

Qxx'i^^'^) / P{x,y]TP)p[y,z;{u- T)f3] 

xp[z,x'; (1 - u)P]V"{y)V"{z) dydz. (54) 

In addition, 2,'(u,t) is continuous in the variables r 
and u and 

QiA^-u,l-r)^Q^^,Ju,r). (55) 

The equation (|49|l remains true for multidimensional 
systems provided that we define the kernel ^,{u, t) as 

in Eq. H4()|l . while for the kernel Q^x'i''^^'^) employ 
the formula 

'3f,x'("''^) = / / dzp{x,y;TP)p[y,z;{u-T)f3] 



xp[z, x'\ (1 — 



^aU^AdlV{y)dlV{z)\). (56) 

i,3 



IV. CONVERGENCE OF THE PA-WFPI 
METHOD 

As argued in Ref. (0), among all possible series repre- 
sentations, the Wiener-Fourier one is special in the sense 
that it has the fastest asymptotic rate of convergence. 
For this reason, it is of interest to establish what the 
asymptotic rates of convergence are for different classes of 
potentials. In this section, we shall show that for poten- 
tials lying in nQ,iy^'2(M), the asymptotic rate of conver- 
gence is always better than ©(l/n^), while for potentials 
lying in nQW^'^(IR), the asymptotic rate of convergence 
is 0(1/71"^). For the latter case, we also establish the 
exact convergence constant. 



A. Potentials with first order Sobolev derivatives 

We begin the convergence analysis with the Kato-class 
potentials that lie in r\aWl^''^{R). For the Wiener-Fourier 
basis, we have 



7„(u,t) = -2 



k=n+l 



smik'Ku) sin(fc7rT) 
P ■ 



(57) 



If f{t) is a square integrable function on the interval [0, 1] , 
then we have the estimate 



ln{u,T)J{u)f{t)dudT 

1 if 

= — ^ — V2 / /(u) sin(fc7rM)d' 

k=n+i L '''^ 
^^r^W E ^ f{u)sm{k 



(58) 



Tru)du 



Now, the first part of Theorem 1 allows us to write 

I3^a^ 



p{x,x';(3)- pI^{x,x';I3) 



< 



E< Xryo{x,x' , (3; a) 



X f duf dTjn{u,T)V'[xr{u)+aBl{d)]{59) 
Jo Jo 



xV'ixrir) +aB"^{d)] 



and using Eq. (|58|l . we obtain 



lim in + If [pix, x'; /?) - pl^{x, x'; /3)] < 



X Ihn E[Xoc{x,x' , P;a)Zn{x,x' , I3;a)] , (60) 

n — >oo 

where 

Zn{x,x\(3;d) 

= V I /" V2sm{k7Tu)V'[xr{u) + aB°{d)]du\ . 
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We notice that the sequence of functions Zn{x, x' , f3; a) is 
monotonicaUy decreasing and convergent to zero. More- 
over, since the functions {^/2sm{k^Tu)}k>l form an or- 
thonormal system, the Bessel inequahty imphes that 

Zq{x, x', f3;a)< f V'[xriu) + aB°{a)]^du. 



We then have the inequahty 
E [Xoc{x, x', (3; a)Zo{x, x' , /?; a)] 

< du dy p{x, y; uP) p[y, x' ; {1 - u)P]V'{y f < oo. 



The last integral is finite because V{x) g HaW^^'^iR). 

In these conditions, the dominated convergence 
theorem^ ^ shows that 

lim E [Xooix, x' , /?; a)Zn(x, x', /3; a)] = 0, 

n — ^oo 

and we have just proved the following theorem 

Theorem 3 // V{x) is a Kato- class potential that lies 
in C\aWa''^{R) and the basis {Xk{u)}k>i is the Wiener- 
Fourier basis, then 

lim{n + lf[pix,x';P)-p^^{x,x';P)] =0. (61) 

n — ^C30 

Observation The previous theorem says that the con- 
vergence is always better than 0{l/n'^) [i.e., it is o(l/n^)], 
but it does not say how much better. If we are interested 
in computing an absolute asymptotic bound, then the 
bound given by the relation 



p{x,x';P)- p'^'^{x,x';(3) 



< 



27r2(n+ 1)2 



xE[Xoo(a;,x',/3;a)Zo(a;,x',;9;a)] < 



27r2mo(n -I- 1) 



(62) 



X d(? / dyp{x,y-ep)V'{yYp[y,x'-{l-e)P] 



is the only choice if no additional information about the 
potential V{x) is available. Theorem 3 remains true as 
stated for multidimensional systems, while in an even 
more compact notation, the estimate given by Eq. H62|l 
becomes 



p{x,x';P)- p^^^{x,x']l3) 



< 



27r2(n-H 1)2 



~e(3H 



E 



{^^Vf 



x' ) d9. (63) 



the class of potentials DaW^'"^ (M.) . Of course, this is also 
true of the Wiener-Fourier series. Moreover, since this 
series is optimal as asymptotic convergence, the rates of 
convergence established in this section set a limit on the 
asymptotic behavior of the partial averaging method. In 
this section, we prove that the asymptotic convergence of 
the PA-WFPI sequence of approximations to the density 
matrix is 0{l/n^) and we establish the corresponding 
convergence constant. 

We start our proof by analyzing the first term in 
Eq. (03. Let us notice that 7„(1 — u, 1 — r) ~ '-fn{u,T) 
because sin[A:7r(l — t)] = — (— 1)*^ sin(fc7ri). By using this 
relation together with Eq. (|39l) . one computes 



f\uf 

Jo Ju 



du / dr gn (w, t)K^ ,^, {u, t) 



1 j-l-U 

du 
Jo 



dr g„ (w, 1 - t)K^^^, (u, 1 - t) 

[ duf dTg„(l-u,l-T)i^f^,(l-M,l-r) 
Jo Jo 



du / dr g„{u,T)Kf^, ^{u,t). 



A little thought shows that the first term in Eq. 149() 
takes the form 



7r2 2^ 



k=n+l 







^, /■" , sin(fc7ru) sin(fc7rr) , , 

du I dr — ^ '-^ -K{u,t), (64) 

l^ 



where 



K{u,t) = kIA^,t)+KP^,Jt,u). 



The asymptotic rate of convergence of this first term is 
dictated by the decay of the terms of the series appearing 
in Eq. H64|) . By integration by parts, we shall later show 
that these terms are of the form Ak/k^ + B^/k'^, where 
{Bk}k>i is a sequence convergent to zero. Because of this 
property, the terms B^/k^ do not contribute to the final 
asymptotic convergence, at least as far as the 0{l/n^) is 
concerned. Indeed, we have 



lim V Bk/k'^ < lim V \Bk\/k^ 



k=n+l 



k=n+l 



< lim sup \ Bk 



k=n+l 



< - hm sup \Bk 

O n^co \k>n 



B. Potentials having second order Sobolev 
derivatives 



As we have shown in Section III.B, all series repre- 
sentations attain their fastest asymptotic convergence on 



which proves our assertion. In the above, we used the 
inequality 

^ n^k^Kn^ / = i 

k=n+l ^ 
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as well as the fact that iJ^ — > implies 



limsup \Bk\ — hm sup \Bk\ I = 0. 

Vfc>n / 

Now, to readily identify which of the various expressions 
make up the terms Bk that decay to zero, the reader 
may employ the Riemann-Lebesgue lemma.^ which in 
our case states that if f{t) is integrable, then 

lim / cos{kTTt)f{t)dt = lim / sm{kTTt) f {t)dt = 0. 

We start with the equation H64|l and by integration by 
parts against the variable r, we obtain 

— / dusin(fc7ru) / sin(fc7rr)_fi'(u, T)dT = 
k Jo Jo 

^ ■ n .K{u,0) ~ cos{kT:u)K{u,u) 
smik-Kuj du (65) 

I r d 

+ — ^ / dusin(fc7ru) / cos{knT)—K{u,T)dT. 
K TT Jo Jo "'^ 

Let us focus on the first term of the relation 1()5|) . By 
integration by parts against u we have 



and 



sm(k7Tu)K{u,0)du 



X(0,0)-(-l)fcX(l,0) 

fc47r2 

d 

cos{kTru)—K{u,0)du (66) 



fc^TT^ Jo ^ ' du 



If . 1 

—r— / sin(A:7ru) cos{knu)K{u, u)du — 

k'^TT Jo 2k-^TT 

X /' siii{2kT:u)K{u,u)du = ^^"'?,7f ^^'^^ (67) 







4fc47r2 



I d 
-. — 7T I cos(2knu)—K(u,u)du, 
'^TT^ Jo du 



respectively. As argued before by means of the Riemann- 
Lebesgue lemma, the last terms in the relations (|66ll and 
(IfTzjl do not contribute to the asymptotic rates of con- 
vergence and therefore, they will be dropped. Moreover, 
even the term K{1, 0)(— l)'^/fc'*7r2 does not contribute to 
the asymptotic rate of convergence because 



E 

k=n+l 



{-If 



k^ 



E 



1 



1 



[n + 2fc + 1)4 {n + 2k + 2f 

1 

X ""dx = - — 

n+l " 



oo „n+2fe+2 

/ 4a;"^da; < 4 / x^^dx 

k=0 ''n+2k+l 



{n + 1)4 



decays faster than Finally, the first term in Eq. H67|) 

is zero because 



K{0, 0) = K{1, 1) = a^p{x, x'- P) [V'{xf + V'{x'f^ 



(68) 



Now, we go back to the formula (|65|l and integrate by 
parts the last term against the variable u. One obtains 

-ir— / dusin(fc7ru) / cos(fc7rr)— if (u, T)dr 
fc-^Ti" Jo Jo OT 

-{-If 



fc47r2 



1 d 1 

cos(fc7rr)— — 1^(1, T)dr + 







di 



d 



X / cos{kTTuf-^K{u,T) 

Jo OT 



du 



r—u 

pa 



k^TT^ 
1 

k^ 



(69) 



X / ditcos(fc7ru) / cos{kiTT)-—-—-K{u,T)dT 

' ' OTOU 



Again, the first and the last terms from the above ex- 
pansion do not contribute to the asymptotic rate of con- 
vergence by the Riemann-Lebesgue lemma. For the last 
term, notice that by symmetry 



pi pu Q2 

/ dMC0s(fc7ru) / cos(fc7rT) K{u, 
Jo Jo o™" 



T)dr 



1 .1 Q2 

cos(knu) cos(kTTT) ^ ^ K(u, T)dwdr 

OTOU 



and then the bidimensional version of the Riemann- 
Lebesgue lemma applies. Therefore, the contribution of 
Eq. H69|) to the asymptotic rate of convergence is due 
solely to the term 



1 



d 

cos( knuf —K(u, r) 

OT 



du 



1 d 

^Jo 



2k ' 

d 

X / cos{2knu)—K(u,T) 
Jo 9t 



2fc47r2 



(70) 



du 



Yet again, the last term of Eq. H7U|) does not contribute 
to the asymptotic rate by the same Riemann-Lebesgue 
lemma. 

Combining everything by means of Eq. (|65|) , we end up 
with the following asymptotic expression for the relation 



37r4n3 

where we also used the equality 



du 



(71) 



lim 71 , 

n^oo \ ^ — ' k 
\ k=n+l 



This last equality can be deduced by letting n — s- cx) in 
Eq. ijBSl of Appendix B. 

The exact expression for the term K{Q, 0) is given by 
Eq. H68|) . To evaluate the term 



du. 
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we need to remember that the density matrix satisfies 
the Bloch equation 

-j^p{x,x';f3)+V{x')p{x,x';l3) = -—p{x,x; (3). 



2mo dx' 



Using the Bloch equation, one can justify by explicit cal- 
culation the following equality: 



-q:^k^^,{u,t) 



2too 



02 

^p{x,y-Tl3) 



92 



xp[y, {u - t)/3] - p{x, y; T(3)-^p[y, z; {u - r) 
xp[z,x'; (1 - u)P]V'{y)V'{z)dydz. 
and by repeated integration by parts against y, 



^K^ (ur) = -^ 



2-^pix,y;TP)V"{y) 
dy 



+p{x,y-Tp)V"'{y) 



p[v,z\ {u-t) 



xp[z,x';{l~u)(3]V'{z)dydz. 

The third order derivative V"'{x) is to be interpreted as 
a distributional derivative. Again in the sense of distri- 
butions, as r ^ u, we have p[y, z; (u — r)/3] S{z — y) 
and therefore, 



as well as the one obtained by permuting x and x' are 
readily obtained by use of the substitution u' — l~u and 
by use of the symmetry of the density matrices. We leave 
it for the reader to utilize these equalities and verify that 



1 d 



du 



mo 



du 

'ji 

(72) 



X / / p{x,y-uf3)p[y,x'-{l-u)P]V"{yfdydu. 



Replacing Eq. H72|) in Eq. (|71|l . one ends up with 
the following expression for the first term of the equa- 
tion (gnj: 

fi2fl3 r ft2^ 
^p{x,x',P)[V'{xf + V'{x'f]+'''' 



2mo 



C (x\e-f''''V"\~^^^-''>"\x') d^j. (73) 



The second term in Eq. H49|l is a little easier to analyze. 
Eq. ljB5|) from Appendix B shows that for the Wiener- 
Fourier series, 



+p{x,y;uP)V"'iy) 



2mo 

p[y,x';{l-u)P]V'{y)dy 



2^p{x,y-uP)V"[y) 
dy 



By integration by parts of the term containing the third 
derivative V"'[y), we obtain 



d 



2mo 



p{x,y;u/3)p[y,x'; (1 - u)/3] 



x^"(2/)'d2;-:^ / V'iy)V"{y) 



■^p{x,y;u(3) 
dy 



^P[y, x'; (1 - u)l3] - p{x, y; u(3)—p[y, x';[l- u)l3\ \dy. 
A similar estimate holds for 



and can by obtained simply by permuting x and x' . The 
equality 



V'{y)V"{y) 



^p{x,y;ul3) 
dy 



xp[y,x'-{l-u)[3]du = j j V'{y)V"{y) 
d 

xp{x',y;ul3)—p[y,x; (1 - u)P]du 
dy 



lim n^gn{u, r) = i^Siu - r) 



in the sense of distributions. Using this equality and the 
relation (|5()|l . we learn that the decay of the last term of 
equation is 



p2 



x''^d0 (74) 



Replacing Eq. and Eq. in Eq. igHl, one ends 
up with the following theorem 

Theorem 4 



lim [p{x,x';P) - pf^'^{x,x';P)] 



p{x,x';p) [V'ixf + V'ix'f 



(75) 



+ / (x e-^«^l/"2e-'^(i-«)^ x') 



>d9. 



The d-dimensional analog can be deduced in a similar 
way by starting with Eqs. (|40|l and (|56|l and by using the 
d-dimensional Bloch equation. It has the expression 
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37r4 

'l27r4 



mo, i Too J 

1,1=1 



) d6'. 



(76) 



We conclude this section by numerically verifying the 
findings of Theorem 4 for the simple case of an harmonic 
oscillator. In fact, we shall verify the following corol- 
lary, which, by trace invariance, is a direct consequence 
of Theorem 4 



Corollary 1 



lim n' 

n — 'oo 



zPA 



(/?) _ f^p{x;P)AZ{x;(3)dx 



Z{(3) 



j^p{x\(})<lx 



where 



^Z[x](3) ^ — — V [x) + 



Stt'^too 



V"{xf 



The Corollary 1 gives an estimate for the relative error 
of the partition function as an average of a suitable es- 
timating function. In practice, such an average can be 
evaluated by Monte Carlo integration if so desired. 

In the remainder of this section, we shall verify the 
Corollary 1 for the simple case of the harmonic oscillator 
V(x) = rriQUj'^x'^ /2. The computations are performed 
in atomic units for a particle of mass too = 1 and for 
the frequency uj — 1. The inverse temperature is set to 
fi = 10. 

Since the density matrix for an harmonic oscillator is 
analytically known,^ the convergence constant for the 
relative error of the partition function can be evaluated 
directly from Corollary 1 to be 



CPA 



ImP{x;/3)AZ{x;P)dx 



= 11.98. 



We can also evaluate it by studying the limit of the se- 
quence 



■Z{fi) 



7PA 



In this respect, we mention that the terms Zf^^{f3) were 
previously evaluated in Appendix B of Ref. (4). Also, 
because the even and the odd sequences and c|^]^ 
have a slightly different asymptotic behavior that gen- 
erates an oscillatory pattern in plots, we choose to plot 
them separately. 

In these conditions, the Corollary 1 says that 



cpA — lim c: 



PA 




FIG. 1: Both the even sequence C2„ and the odd sequence 



C2n+i are seen to converge to cpa 
predicted by Corollary 1. 



11. 



which is the value 



and the prediction is well verified for the harmonic oscil- 
lator, as Fig n shows. We regard this numerical example 
as strong evidence that the analysis performed in the 
present paper is correct. Indeed, it is hard to believe 
that the exact agreement between theory and numerical 
analysis for the harmonic oscillator is accidental. 



V. SUMMARY AND DISCUSSION 

In this paper, we have performed a complete analy- 
sis of the asymptotic rates of convergence of the partial 
averaging sequence of approximations for the density ma- 
trix. We have found that there are two natural classes 
of potentials on which different series may achieve their 
fastest convergence. If the two-point correlation function 
7„(m, r) for the tail series is positive, then the natural 
class is r\aW^''^{M.) as shown by Theorem 1. More gen- 
erally, all series representations achieve their fastest con- 
vergence on the space of potentials r\aW^''^{M.). This is 
the statement of Theorem 2. Moreover, in both cases, 
the cited theorems provide the exact asymptotic rates of 
convergence for arbitrary series. 

In Section IV, we analyzed the special case of the 
Wiener-Fourier series representation. Theorem 3 asserts 
that the asymptotic convergence for this series is o{l/n'^) 
on the Kato-class potentials that lie in r)aW^''^{R). 
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Moreover, the equations H62|l and (|63|l provide useful 
(though conservative) bounds that control the absolute 
asymptotic error. The Wiener- Fourier series achieves 
its fastest asymptotic convergence on the Kato-class po- 
tentials that lie in nQ.W^'^(R). This convergence is 
0(l/n^), with a convergence constant given by Eq. H75(l 
for monodimensional systems and by Eq. H76|) for multi- 
dimensional systems, respectively. 

The present results provide a relatively complete char- 
acterization of the asymptotic convergence characteris- 
tics of the partial averaging methods. Beyond establish- 
ing sharp estimates of the convergence constants for ap- 
plications where the potential energy has second-order 
derivatives, it should be possible to utilize the present 
methods to characterize the asymptotic behavior of other 
properties. In particular, we would speculate that it 
should be possible to obtain convergence constants for 
both the H-method and T-method energy estimators us- 
ing arguments similar to those developed in the present 
paper. Moreover, we expect that the theorems proved 
in this paper will play a crucial role in establishing 
the asymptotic rates of convergence for the reweighted 
method, a technique related to the partial averaging 
method. 
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APPENDIX A 

In this appendix, we compute E(i?")^, E(_B^)^, and 
E(i?°i3^). We begin with the last expression. We have 



Jo 



I - ^Afc(u)Afc(T) 

k=l 

\k{u')\k{T')dT'du' 



(Al) 



fe=l 



/" r [5{u' - t') - 1] dr'du' 
Jo Jo 



where we used the completeness relation. In this sense, 
the reader should remember that {Xk{u)}k>i together 
with the constant function make up a complete basis in 
L^([0, 1]). Now, let Iu{u') and It{t') denote the indicator 
functions of the intervals [0,m] and [0, r], respectively. 
The last integral in Eq. (|A1|) becomes 

E(B°B^) = /" [ Iu{u')Ir{T')S{u' -T')dT'du' -UT 

Jo Jo 

= / Iu{u')lT{u')du' — UT — min{u,T) — UT. (A2) 
Jo 



Then, the expressions for E(i?")^ and E(S")^ are read- 
ily obtained by setting w = r in (|A2I) . We notice that 
the above expectations are independent of the particular 
random series representation of the Brownian bridge, as 
they should be. 



APPENDIX B 

In this appendix, we prove that for the Wiener-Fourier 
series we have 

hm Io\^f^rnMdudr ^ /\(,,,)d. (Bl) 
/o /o 7«(",T)2dwdT Jo 

for all continuous functions h{u, r). This is the statement 
of Eq. I|13|l in Section II. B. We start with the equations 



sin(fc7ru) sin(fc7rr) 



(B2) 



k=n+l 



and 



pi ri 1 °° 1 

Jo Jo 

By letting n — s- cx) in the sequence of inequalities 

1 roc oc 

"'"+1 k=n+l 



3 (n+ 1)3 



one deduces that 



< I x^'^dx = i, (B3) 



1 rl 



lim n3 / / ^^(u,T)2dudr = -i-. (B4) 
Jo Jo Stt'' 

Therefore, Eq. (|Bip can be reformulated as 

»1 rl rl 



Sn^ lim 



"'0 



7„(w, r) h{u,T)dudT = j h(u,u)du 

(B5) 



We prove the relation (|B5|) in two steps. The first step 
is to set h'{u,T) — t) — h(u^u)\ and show that 

3tt^ \im [ [ j„{u,T)^h'{u,T)dudT ^ 0. (B6) 
Jo Jo 

For that purpose, pick an arbitrary r; > 0. Let us notice 
that h'{u,T) is continuous on the compact set [0,1] x 
[0, 1], thus bounded by a constant M < oo and uniformly 
continuous. The second property implies that there is 
1 > e > such that |/i'(u',r') — h'(u,T)\ < -q whenever 
(m'--w)^-I-(t'— r)^ < e^. Now, notice that by construction 
h'{u,u) = for all < u < 1. Since any point (w, r) in 
the set 

= {(u,t) G [0,1] X [0,1] : |w-r| < e} 
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is contained in the ball of radius e centered about 
the point (u, u), by uniform continuity, it follows that 
h'{u, r) < 77 on I^. 

Let us break the integral in Eq. (|B6|) in two parts: 
one over the set and the other on its complementary 

— [0, 1] X [0, 1] — le- For the integral over the set 1^ we 
have 



Stt** lim sup 71'^ / / jniu,T)'^h'{u,T)dudT 

n— ^00 J Jj^ 



< 



Stt^ lim sup Ti"^ / / 7„(u, rj'^dudr 



< V- (B7) 



For the integral over the set I^, we have the estimate 
Stt* lim sup n'^ / / 7„(m, r)^ft,'(u, T)dudr 

)l— ►CXD J Jjc 



< M 



Stt limsupn / / 7„(u,t) dudr 

n^oo J J Jo 



(B8) 



We intend to show that the last limit in Eq. ljB8|l is zero 
and for this purpose, we need to construct a sharp bound 
of the function 7„(u,t) on the set I^. It turns out that 
we first need to study the behavior of the functions 



E 

k=n+l 



cos{kTru) 
fc2 



on the interval (0, 2) 
We have 



. / /„^ cos(fc7ru) 

k=n+l 

1 ^ sin[(fc - l/2)7rii] 



5 Y. 



sin[(fc + l/2)7rw] 



fc2 



k=n+l 



k=n+l 



k=n+l 

1 sin[(n + l/2)7rw] 
'2 (n + l)2 

sin[(fc + l/2)7ru]. 



(B9) 



From Eq. HB9|) we learn that 



sin(7r?i/2) 

00 



k=n+l 



cos(fc7ru) 



fc2 



< 



1 



fc=n+l 



1 1 



fc2 (fc + 1)2 



2(n+ 1)2 
1 

^ (n + l)2- 



Therefore, 



cos(fc7ru) 



^ fc2 



< 



where the function 



F{u) 



sin(7ru/2) 



(BIO) 



;F{U), (Bll) 



is continuous on the interval (0, 2), decreasing on the in- 
terval (0, 1] and increasing on the interval [1,2). 

Using the bound given by Eq. 1B11|I together with the 
inequality (a — 5)2 < 2(a2 + Ip'), one computes 



< 



OC 

E 



cos[A:7r(w — r)] — cos[A:7r(w + t)] 



/c2 



cos(fc7r|'u — t|) 



,k=n+l 



E 

.k=n+l 



cos(fc7r|u + t|) 

fc2 



(B12) 



< 



^^AF{\u~r\f+F{\u + r\f] 



On the set we have e < |m — t| < 1 and e < |w + t| < Corroborating with Eq. 1B12|I . we obtain 
2-e. Therefore, F{\u-t\) < F{e) and F{\u + t\) < F{e). ^ 

ln{u,T)' < F{ef if{u,r)el^. (B13) 
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Replacing the estimate (|B13|I in Eq. (jB8|l . one readily 
concludes that the limit in Eq. (|lj8|l is 0. Adding the 
limits of Eq. I)B7|I and Eq. (|B8|I , we conclude that 



37r'*limsupn^ / / 7„(u, t)^/i'(u, T)dudT < 77. (B14) 

n— ►oo Jo Jo 

Since 77 > is arbitrary, the equality (|B6|) is demon- 
strated. 

In the second step of the proof, we show that Eq. ljB6|l 
implies the relation (|B1|) . From Eq. (|B6|) . one deduces 
that 

Stt^ lim / / 7„(m, t)^/i(w, r)dMdr 
Jo Jo 

= 3tt^ lim n^f [ jn{u,Tfh{u,u)dudT. (B15) 
Jo Jo 

By explicitly computing the integral over t, the last re- 
lation becomes 



Stt* lim / f„{u)h{u,u)du, (B16) 

n — ^00 



where 



2 ^ sin(fc7ru)^ 1 

k=n+l 



- F 



1 — cos(2fc7ru) 



A;=Ti+l 



By continuity, h{u, u) is integrable and if we set 



Bk = I cos{2kTTu)h{u,u)du, 
Jo 



(B17) 



the Riemann-Lebesgue lemma says that 

lim \Bk\ = lim sup \Bk\ — 0. 

We notice that 



Stt lim 

n — >oo 7r4 



^ /"^ cos(2fc7ru) 



< 37r'' hm 



fc=n+l 



00 \ 

^ sup \Bk\ 



■ A;— 71+1 / — 

< lim sup \Bk\ = 0. 



Therefore, the terms containing cos(2A;7ru) in the 
Eq. ljB17p may be dropped and the quantity ljB16|l be- 
comes 



\k=n+l / ■^'^ 



u — I h{u, u)du, 


(B18) 

where we used (|B4(l to compute the last limit. The rela- 
tion (|B1|I is therefore proved. 

APPENDIX C 

It is well known that if A, i? > 0, and a = Cj^ AB 
such that \a\ < 1, then the following Mehler's formulaiS4 
holds for all / and g whose squares have finite Gaussian 
transforms 



[/f;W(xo,yo) = y^dxy^d7;-^_=exp(^-- ^^-^^ j /(^o + x^j/o + 2/) 

= / dx\ dy^^=l==expf-l "^' +/^'~f^" )/(a:o + x^/^)g(7/o + ^/v^) (CI) 
27r VI - a2 V 2 1 - / 

= ^E"' f dx f dye-^''"+y"^/^Hk{x)Hk{y)f{xo + xVA)g{yo + yVB). 



In the above, the functions Hk (x) are the normalized Her- 
mite polynomials corresponding to the Gaussian weight 



dfi{x) — 




They form a complete orthonormal basis in the Hilbert 
space L^(]R), which is endowed with the scalar product 

im = / i;{x)cj,{x)dfi{x). 
Jr 



Let us notice that according to our hypothesis, the func- 
tions /(xo + x^/A) and g{yo + x\/B) as functions of x are 
square integrable against d^{x) and thus they lie in the 
Hilbert space 

By repeated integration by parts, the formula (|C1|) is 
shown to equal 



[f9]ABci^o,yo) - E -]^fAM/B\yo), (C2) 

fc=0 
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where in general (xq) is the fc-order derivative of 
1 



3-^ /(2^)/(xo + z)dz. 



Let us notice that the series (|C1|) can be extended to the 
case a = 1, too. Indeed, the last series in Eq. (|Cip for 
the case a = 1 is nothing else but the Bessel series 

OG 

k=0 

which is convergent to 

(/(xo + -VI)|.9(yo + -/B)) 

= / f{xo+xVA)g{yo + xVB)dii{x). 

Next, we proceed to establish the relation p8|l . We 
start with the identity 

En [Uoo{x,x' , P; d) - En C/oo(a;,x',/3;a)]^ 

= E„ U^{x, x', (3- af - [E„ U^^, x\ /?; a)f . (C3) 

Using the notation introduced in Section II. A, we have 

E„;7oo(a;,x',/3;a) =F„,„[a;,(u)+a5;'(a)]. (C4) 

Moreover, 



"^nUoaix, a;', /?; a) 



1 .1 

du I dTE„y[xr(it) 



+aS:,{d) + aB-^{d)]V[xr{T) + aS^:{d) + aB^:{d)] (C5) 

and remembering that the variables B"{d) and B"{d) 
have a joint Gaussian distribution of covariances given 
by Eq. O, we obtain 



EnUoo{x,x',(3;a)'^ = f du j dr ( dx [ dy ^ 

Jo Jo Jm Jr '^T^^n\ 



where V ^ ^(x) is the k order derivative of Vu,n(x). With 
the help of Eq. (jC4|l , one recognizes the first term of the 
above series to be [E„[/oo(a;, a;', /3; a)]^, so that Eq. HC3|I 
becomes 



En \Uoo{x,x' a) - E„C/oo(a;,a;',/3;a)]^ 

^Y.u dTG„(zi,T)'= (C6) 

fe=i Jo Jo 



As mentioned in Section II. B, the kernels 7„(u, t)'^ and 
therefore G'„(u,t)'^ are positive definite for all integers 
fc > 1 so that the terms of the series (|C6|I are all positive 



I.e. 



\u C dTGn{u,TfV^^l[Xr{u) + CjSl{d)] 

Jo 

V^^l[xr{T) + aS^^{d)]>Q 



for all fc > 1. Because of this property, the series (|C6|) 
will plays the fundamental role in establishing the various 
results in this appendix. Truncating the series llCfill to 
the first term, we obtain the inequality 



\Uqo{x,x ,(3] a) - E„ Uoo{x,x' , (3; a)f > E„T,' (a;, x', /5; a), 



(C7) 



where 



(m,t) 



X exp 



l x^rl{T)+y^rl{u)-2xyGniu,T) 

xV[Xriu) + (75^(5) + x]V[xriT) + aSl\d) + y] 
where Gniu^r) — a'^^n{u,T) and 

Al{u,T)^Tl{u)Tl{T)~Gn{u,Tf. 



Using the expansion (|C2|I . one may write the above 
integral as the sum of the series 



,x',/?;a)2 = f^l f du f dTGn{u,T)^ 
k=o Jo Jo 



xvi%riu)+aS:i-a)]V^'l[xr{T)+aS:^id)l 



Tn{x,x',(3;a) ^ du drG„(u,T) 
Jo Jo 

xvi%r{u) + aS:{-a)]vi'i[xr{T) + aS^id)]. 



In a similar fashion, if we truncate the series ifUe)! to the 
first two terms, we obtain 



En [Uooix, x', P; a) - E„ Uooix, x' , /?; a)] > E„y,'(a;, x' , l3; a), 

(C8) 

where 



Yn{x,x',(3;a) = Tn{x,x' , /3;a) + i / / dTGn{u,Tf 

^ Jo Jo 



Now, we assume that V{x) S HaW^'^iR). The in- 
equality 1/kl < l/{k— 1)1 and the positivity of the terms 
of the series (|C6|I imply 
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E,,[Uooix,x',(3;a)-K„Uo.ix,x',P;af < / du dTG„(u,T) V-G„(w,t)'= 

xVi't'\xr{u) + aS:{a)]vi'^'\xr{r)+aS-{a)]'j = j\u j^rG^M 
xE„ { [xr{u) + (TBO(a)]y(i) [xr{T) + aB°(a)]} = E„r„(:E, x', /3; a), 



(C9) 



where 



where 



1 .1 

du/ ATGn{u,T)V^'^^[xr{u) + aBl{a)] 
Jo 



xy(i)[x,(T) + ^S°(a)]. 



Similarly, by employing the incquahty 1/fc! < l/[2(fc— 1)!] 
if fc > 2, one proves the sharper inequality 

En Poo{x,x\ P; a) - E„ Uooix, x' , (3; a)f < E„y,"(a;, x' , f3] a), 

(CIO) 



Y;:{x,x\P; a) = -[T:^{x,x',f3;a) + T,,{x,x',l3;a)] 



If in addition V{x) e HaW^^^C 



K(zi) + a5,7(a)]y^+'^K(r)+a5,"(a)] < - / du drG„(w,rf K]-G„(^i,r)* 



E, 

(fe+2) 



xv['t'\xr{u) + aS:{a)]V^'^^\xr{T) + aS^ia)] 



r 



(Cll) 



One recognizes the last series in Eq. (jClljl to be equal to 
E„ [V^^^ixriu) + aB°ia)]V^^^[xriT) + <jB°ia)]} 



The multidimensional case is left to the reader and 
can be obtained starting with the following d-dimensional 
analog of Eq. ljC6|l . Let Ik be the set of all vec- 
tors k = {ki, . . . , kd) of positive components such that 
-Vu!n[xr{u) + aSl^{a)]vl^l[xr{T) + a5;?(a)] (C12) /ci + • • • fed = fc and let ... ^^J{x) denote the respec- 

tive partial derivative of the d-dimensional function f{x). 

and therefore, the inequality (|45|l of Section III.B is Then, 
proved. 



— 1 

k=i ' •'^ "^o 



E 



kl 



ki\---kd\ 



• --^T {dL...,k,VuA^r{u)+aS-{a)] dl_k,V rA^^r) + aS-{a)]] \ . (C13) 



The equation ljC13|l can be deduced by remembering that pendent processes. Then, one successively applies the 
the Brownian bridges B^ f. for each dimension are inde- 
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expansion (|C6|I for each dimensfon to obtain Eq. ljC13|) . 
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